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Abstract 

We give an adaptive, fixed-point version of Grover's algorithm. By this we mean that 
our algorithm performs an infinite sequence of gradually diminishing steps (so we say 
it's adaptive) that drives the starting state to the target state with absolute certainty 
(so we say it's a fixed-point algorithm). Our algorithm is motivated by Bloch sphere 
geometry. We include with the ArXiv distribution of this paper some simple software 
(Octave/Matlab m- files) that implements, tests and illustrates some of the results of 
this paper. 
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1 Introduction 



Grover proposed his original algorithm in Ref. [T]. His algorithm takes a starting state 
towards a target state by performing a sequence of equal steps. By this we mean that 
each step is a rotation about the same fixed axis and by the same small angle. Because 
each step is by the same angle, the algorithm overshoots past the target state once it 
reaches it. Grover later proposed in Ref.|i2] a "vr/S fixed-point" algorithm which uses 
a recursion relation to define an infinite sequence of gradually diminishing steps that 
drives the starting state to the target state with absolute certainty. 

Other workers have pursued what they refer to as a phase matching approach 
to Grover's algorithm. Ref.fS] by Toyama et al. is a recent contribution to that 
approach, and has a very complete review of previous related contributions. 

In this paper, we describe what we call an Adaptive, Fixed-point, Grover's 
Algorithm (AFGA, like Afgha-nistan, but without the h). Our AFGA is motivated 
by Bloch sphere geometry. Our AFGA resembles the original Grover's algorithm in 
that it applies a sequence of rotations about the same fixed axis, but differs from it in 
that the angle of successive rotations is different. Thus, unlike the original Grover's 
algorithm, our AFGA performs a sequence of unequal steps. Our AFGA resembles 
Grover's 7r/3 algorithm in that it is a fixed-point algorithm that converges to the 
target, but it differs from the 7r/3 algorithm in its choice of sequence of unequal 
steps. Our AFGA resembles the phase-matching approach of Toyama et al., but their 
algorithm uses only a finite number of distinct "phases" , whereas our AFGA uses an 
infinite number. The Toyama et al. algorithm is not guaranteed to converge to the 
target (in the single-target case, which is what concerns us in this paper), so, it is not 
a true fixed-point algorithm. 

This paper was born as an attempt to fill a gap in a previous paper, Ref.[l], 
which proposes a quantum Gibbs sampling algorithm. The algorithm of Ref. |1] re- 
quires a version of Grover's algorithm that works even if there is a large overlap 
between the starting state and the target state. The original Grover's algorithm only 
works properly if that overlap is very small. This paper gives a version of Grover's 
algorithm without the small overlap limitation. 

We include with the ArXiv distribution of this paper some simple software 
(Octave/ Mat lab m-files) that implements, tests and illustrates some of the results of 
this paper. (Octave is a freeware partial clone of Matlab. Octave m-files should run 
in the Matlab environment with zero or few modifications.) 

Some vocabulary and notational conventions that will be used in this paper: 
Let Ns = 2^^ be the number of states for Nb bits. Unit vectors will be denoted 
by letters with a caret over them. For instance, f. We will often abbreviate cos(^) 
and sin(6'), for some 6* G M, by Cq and Sg. 

We will say a problem can be solved with polynomial efficiency, or p- 
efRciently for short, if its solution can be obtained in a time polynomial in Nb- 
Here Nb is the number of bits required to encode the input for the algorithm that 
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solves the problem. 

By compiling a unitary matrix, we mean decomposing it into a SEO 
(Sequence of Elementary Operations). Elementary operations are one or two qubit 
operations such as single-qubit rotations and CNOTs. Compilations can be either 
exact, or approximate (within a certain precision). 

We will say a unitary operator U acting on can be compiled with 
polynomial efficiently, or p-compiled for short, if U can be expressed, either 
approximately or exactly, as a SEO whose number of elementary operations (the 
SEOs length) is polynomial in Nb- 



2 Review of Pauli Matrix Algebra 

The Pauli matrices are 
ax = 

Define 
For any a G M^, let 



' 


1 " 




' -i' 




" 1 " 


1 





, cry = 


i 




-1 



o" = (o"x, cry, 



= a ■ a . 

(I often refer to "Paulion", related to a quaternion). 

Suppose a, 6 G M^. The well know property of Pauli Matrices 

(^p^^q ^P,q i^pqr^^r i 

where indices p, g, r range over {1, 2, 3}, immediately implies that 



(1) 

(2) 
(3) 

(4) 
(5) 



A A ^A ,A A 

K r = r-2a{a-r) 




Figure 1: Kaf is the reflection of f on the plane perpendicular to a. 
For any unit vectors f, a, define 



= r — 2a{a ■ f) . 



(6) 
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From FigJH we see that Ka is a reflection: it reflects the vector r on the plane 
perpendicular to a. Furthermore, —K^ is a pi rotation: it rotates the vector f by an 
angle of pi, about the axis a. Note that 



(a ■ f + zcraxf)cra 



(7) 



\a-r\a. 



a 



(axr)xa 



{a ■r)aa- (yr^a{a-r) 



(9) 
(10) 




Figure 2: A pi rotation of z about dr yields f. 



The eigenvectors of az are given by 



az\±z) = (±1)1 



where 



1 




|o), 1-^) 





1 



fill 



:i2) 



It is also convenient to consider the eigenvectors of af, for any unit vector r. Given 
f, we can always flnd a unit vector dr such that f is obtained by pi rotating z 
about dr- (See FigJ2]). Eg. ( ITT]) implies that 



(o-a^o-iO^aJo-ajO) = (Ta,|0) . 



Furthermore, 



Hence, if we deflne |r) by 



(13) 



(14) 



\r) = aar\0) 



(15) 
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the Eq. (fT3!) becomes 



af\f) = \r) 



Eq. f|T6l) itself implies that 



(16) 



(17) 



(Jrl ± f) = ±1 ± f) . 

Thus, I ± f) defined by Eq. (fT5l) are the eigenvectors of af with eigenvalues ±1. 

Given f in polar coordinates, we can express |f) in terms of the polar coordi- 
nates of f, as follows. Suppose 



sin 6 cos I 
sin 6 sin < 

cos 6' 



where 6 G [0, |] and G [0,27r). Then, by definition, 
9 coordinate: 



has the same 6 but half the 



sm 



cos< 



smi^^J sm( 
cos(^) 



Thus, 



(19) 



l^) = ^a.|0) 
From Eq.f l20|) . we get the following table. 







cos(f) 


Qiijt 1 lj(Xy 




e^* sin(|) 



(20) 









\±x) 


1 

~ V2 


1 
±1 






l±y) 


1 

~ V2 


1 

±i 





1^) 



1 





(21) 



In general, the assignment f i— |f) yields a map of into C^, the 2-dimensional 
complex vector space spanned by complex linear combinations of |0) and |1). If we 
confine ourselves to the x — z plane of M^, then that plane is mapped into the half- 
plane of all real linear combinations of |0) and |1) with non-negative |0) component. 
(Kets that differ by a phase factor or a normalization constant are equivalent). 

Just like it is useful to consider the projection operators |0)(0| and |1)(1|, it is 
also useful to consider the projection operator |r)(f|. Recall the usual definitions of 
the number operator n and its complement n: 

n = Pi = |1)(1| = , (22a) 
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n=l n = Po = |0)(0| = ^ + ^^ . 


(22b) 


Rotating the coordinate system so that z goes to f, we get 




= = 1 r){ r\= 2 


(23a) 


1 + (T- 

rir = 1 — JT-r = -Pr = = ■ 


(23b) 


Note that if we define the refiection operator K\f^ by 




X|,) = l-2|f)(f|, 


(24) 


then 




fsTl^) = 1 - Irif = (-1)"^ = -af . 


(25) 


In general, the vectors in arc "packed twice as densely" as the corresponding 
vectors in M^. By this we mean that the angle between two vectors in is always 
half the angle between the corresponding vectors in R^. To prove this, note that 


I(ri|r2)|^ = tr(|fi)(fi||f2)(f2|) 


(26) 


= |tr[(l + (j,J(l + (7,J] 


(27) 


1 , 


(28) 


= ^tr(l+ri-f2) 


(29) 


= -(l+ri-r2). 


(30) 


Thus, if fi ■ f2 = cos a, then (fi f2) = \J^^^^^ = cos(|) . 

Suppose ^ e R. A Taylor expansion easily establishes that 




e'^'^^ ^ cose + iaz sine . 


(31) 


Rotating the coordinate system so that z goes to f, we get 






(32) 



Suppose that a and b are two unit IR vectors which make an angle e between 
them. Let h — be the unit vector normal to the plane defined by a and b. 

\axb\ ^ 

Then 
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Figure 3: A finite rotation of f about axis a by an angle of ^. 



CaO"b = a ■ b + icr^^i^ = cos 6 + iaa sin 6 = e^^"" . (33) 

Thus, any SU(2) element e^^"", where 6' G M and n is a unit vector, can be 
expressed (non-uniquely) as a product of two Paulions aa and cr^. 

From FigJ31 it is clear that if Rd{C) is the rotation operator that rotates any 
unit vector f, by an angle ^ G M, about an axis defined by the unit vector a, 
then 

RdiO^ = ^(^ ■ ^) + sin(^)a X f + cos(,^)[f — a(a • f)] . (34) 




Figure 4: An infinitesimal rotation of f about axis a by an angle of A,^. 

If f, a are unit vectors and is an infinitesimal real number, then 

e 2 '^afc 2 ° ^ - z— ((TaCTf - cr^^aa) = crr+AC(axf) • (35) 

From FigJU f + A^(a x f) is just the vector f after an infinitesimal rotation i?a(A^). 
By applying successively a large number of infinitesimal rotations Ra{^^), we get a 
rotation -Ra(0 o'^^^ ^ finite angle ^ G M: 

e-4-^a,e^i'^^ = . (36) 
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3 A SEO That Takes |s') to \t) 

In this section, we will explain a formalism and accompanying geometrical picture 
that can be used to describe both the original Grover's algorithm, and the AFGA 
proposed in this paper. 

Our goal is to find a SEO of SU(2) transformations that takes a starting state 
\s') e to a target state \t) e C^. Without loss of generality, we will take 



1^) 



|0) 



and 



|s') = \s') where s' — {sm^)x + {cos^)z 



(37) 



(38) 



with 7 e [0, tt]. 

Let Ej for j = 0, 1, . . . , jmax denote 



(39) 



where e R and the aj are unit vectors. Suppose is generated as follows 

j—jmax ^0 

(The arrow in the subscript of the product sign indicates the order in which to multiply 
the terms, this being an ordered product.) It follows that 



(oi( n ^^oi^') 

j —jmax ^ 



(oi( n n 4)io) 

j—jmax ^0 J— ^jmax 

(oi( n ^^■)(^)( n 4)io) 

j=jmax-^0 ^ ^ o=n-^o 



(41) 
(42) 



j — yjmax 



^(i+(oi( n ^"'""'^M n e^^-.)io)|43) 

\ j—jmax ^0 j—0 yjmax / 



(44) 

(45) 



If we define ERR by 



ERR = 1 



(oi( n E,)\s') 

j—jmax ^0 



(46) 
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then Eg. (14 5 p imphes that 



ERR=-{l-[sf,. 



(47) 



ERR > is a measure of error. By Eg. (1461) . ERR decreases towards zero as the SEO 
takes |s') closer to \t). Eg. ( 1471) agrees with our expectation that ERR goes to zero as 
the Z component of Sfm approaches one. 

In the original Grover's algorithm, |s') = U\s)^ where \s) = |0)®^^ and U = 
jj(^Nb ^ (Here H denotes the one-bit Hadamard matrix -^{crx + <^z)- Also \t) = 

\t]\fg_i) (g) . . . \ti) \to) with ti = 61° for some Zq G {0, 1, . . . , A''^ — 1}. Furthermore, 



n 



QNstepa 



(4^ 



j — 3 max ^0 



where N steps is the number of steps (i.e., gueries) and 



G 



_(_1)K>(^'I(_1) 



\t){t\ 



(49) 
(50) 
(51) 



Note that G corresponds to a small rotation about the y axis. Indeed, 

G = -(-l)^^'(-l)"^ 

= —as'^z 

= -{§' ■ z + ias'xz) 

= —(cos 7 — icTy sin 7) 

— g«('f-7)o"y — g*^'^!' 

where A7 = 2(71 — 7) ^0- Thus 

Rs'{-7T)R,{-n) = R^{-A-f) . 



(52) 
(53) 
(54) 
(55) 
(56) 



(57) 



From cos(| 



{s'\t) = it follows that A7 = 0{^). 
Note also that G can be p-compiled trivially. Indeed, if Pq and Pi are defined 
as in Egs.(l22]), 



(-i)i*)(*i = (_i)n^fo"'^*.(/3) ^ 



(58a) 



and 



JO^s)(0^s| 



jj(^Nb _ jj<^Nb 



-l)n^So' PoW)] H^Nb _ (58b) 
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Next, let's describe our AFGA. We take 



j=jmax^O j=OO-5>0 

= Yl (e5("^+^^)e*^"=-'e*^"^) , (60) 
j =00—^0 

for some AA G [0, vr] and some infinite sequence of aj G M. In the original Grover's 
algorithm, A A and all the aj are fixed at vr. That's why we say our algorithm is 
adaptive. We will also assume that 

Rs'{-aj)R,{-AX) = Ry{-Ajj) (61) 

for some A'jj G M, and that we know how to p-compile e*"j''*'^^'*'' and e'^'**'*^^*'. If 
\s') and |t) are the same as in the original Grover's algorithm, then we do indeed 
know how to p-compile these operators. (Just replace the (—1) = e^^ phase factors 
in Eqs.dSS} by e*°^ and e'^^.) 

Our AFGA has a total of two input parameters: AA G [0, vr] (the angle of each 
consecutive z rotation), and 7 G [0, vr] (the angle which s' makes with z). For the 
AFGA to be fully specified, we still need to specify, as a function of these two input 
parameters, a suitable sequence of aj that makes ERR converge to zero. We will do 
this in the next section. 



4 Bouncing Between Two Longitudes 

In this section, we give a suitable sequence {aj}°^Q for our AFGA, as a function of the 
two input parameters AA and 7. We will do this guided by the geometrical picture 
of bouncing between two great circles of longitude of the unit sphere. 

It is convenient to define two sequences of unit vectors {sj}'jLQ, and {rj}°^o 
by the recursion relation: 

So = s' (62) 
rj = R^{-A\)sj, 

Sj+i = Rs'{-aj)rj (63) 
%+i = Ry{-A'yj)sj 

for J = 0, 1, . . .. According to Eqs. fl63l) . 



{fj,Sj) = AA, 



Z{sj+i,fj)=aj, . (64) 
Z{sj+i,Sj) = A7j 
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A 

-z 



A 

-z 



Figure 5: Unit vectors and angles used in our AFGA. Note how we "bounce 
between two longitudes" . 



See Figsj5]for a geometrical picture of the Sj and fj sequences of vectors. The 
arrowheads of all the Sj vectors lie on the great circle of longitude located at the 
intersection of the x — z plane and the unit sphere. The arrowheads of all the fj 
vectors lie on the great circle of longitude which makes an angle of AA with the Sj 
great circle of longitude. 

The angles 7^ and hr^j are defined in terms of the vectors Sj as follows: 



A 

s. 



j+1 



A 

r. 




Figure 6: Geometry defining angle aj and radius L. 




(65) 



and 



Z(sj, z) = -fj , Z{sj+i, z) = 7j+i A7, 



A7j- = 7j - 7j-+i 



(66) 



for j = 0, 1, . . .. 
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Note from Fig|5] and |6] that 

Z(fj, s') = Z(%+i, s') = 7 - Ij+i = 7 - 7i + ^Ij 



Since 



fj = Rz{—AX)sj 

= z{z ■ Sj) - z X sjSax + [sj - {sj ■ z)z]Cax 
= zC^^{l - Cax) - vS^^Sax + sjCax , 



(67) 



(68) 
(69) 
(70) 



it follows that 



fj ■ s' = CjC^^{l - Cax) + cos(7 - ^j)CAx 
= C^C^. + S^S^.C AX 



(71) 
(72) 



As is usual in the C programming language, define 9 = atan2(?/, x) only if 
tan 6' = y/x. (Think of the comma in atan2(?/, clS db slash indicating division). 
Eq. fl671) implies that 



7 — 7j + A'jj = atan2(±-v / 1 



{rj ■ s'y,rj ■ s') 



(73) 



Assume /L{fj,s') G [0,7r]. This means we choose the solution with the positive sign 
in Eq.(I73l). 



To summarize, we've shown that 



where 




fj'S — C^C^. + S^S^.C AX 



(74b) 



oo 
jSj=0 



Eqs.( !74|) and 7o = 7 allows us to find the sequence {7j}jLo- The sequence {7^} 
and So = s' allows us to find the sequence {sj}^^^. 

Next we solve for the sequence {aj}°2,g, assuming that we know the two input 
parameters AA and 7, and the sequence {7j}°lo- ^^"^ ^3 ^"7 considering FigEl 

That figure defines L as 



L = \rj-s\s' ■rj) \ = ^1 - (r^- ■ s'f = \ sin(7-7j- + A7j)| = sin(7-7j- + A7j) . (75) 
Note also that 



It follows from Figj6]that 



(76) 
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sin aj = [rj - s\s' ■ rj)] x [sj+i - s'{s' ■ rj)] ■ {-§') (77) 
= fjXSj+i-{-s') (78) 
= '^A7, [rj X X • s] + Ca^^ [-r^ x Sj ■ s'] . (79) 

But by virtue of Eq. (!70|) for fj, 

fj X {y X Sj) ■ s' = -{sj ■ s){fj ■ y) = cos(7 - 'yj)S^^SAx , (80) 

and 

f-j X Sj ■ s' = -{y X Sj ■ s')Sy^SAx = - sin(7 - 7j)5'^,5'aa (81) 

so 

sinoj = sin(7 - 7^ + A7j)S'^^.S'aa • (82) 
It also follows from FigEJtliat 



cos aj = [rj - s{s ■ Vj)] ■ [sj+i - s{s ■ Vj)] 
= [rj - s'{s' ■ fj)] ■ Sj+i 
= [fj - s{s' ■ fj)] ■ [-y X %^A7, + sjCa^.] 

^A7, [^i X % ■ y + sin(7 - -fj)fj ■ s'] 
+Ca7, [rj ■ Sj - cos(7 - -fj)fj ■ s'] 



^3) 
M) 
55) 

^6) 



But by virtue of Eq. fITU]) for f 



J' 



fj X Sj ■ y = {z X Sj ■ y)C^^{l - Cax) = Sj^C^^{l - Cax) 



and 



(87) 



frs,=C'(l-C, 



AX 



Cax — + S'^.Cax 



so 



2 / '^A7, [C-ij 5'7, ( 1 - Caa) + sin(7 
" I +Ca7, [CI + ^^,Caa - cos(7 - 


- 7i)^j ■ s'] 

- 7i)^i ■ s'] 


(89) 


To summarize, 


aj = atan2(L^ sin aj,L'^ cos aj) 




(90) 



where sin aj and cos aj are given by Eqs. (J82l) and (IHS 
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5 Software and Numerical Results 

In this section, we describe some simple software that calculates, among other things, 
the phases aj used in our AFGA. We also present and discuss some examples of the 
output of the software. 
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1.4921e-03 


1.4921e-03 


1. 


OOOOe+00 


-2. llOle-03 


-1.6948e-15 


1 


.OOOOe+00 


20 


8.6014e-02 


5 


. 1447e-01 


-1 


.0615e-03 


-1 


.0615e-03 


1. 


OOOOe+00 


1. 


5012e-03 


1.1601e-16 


1 


.OOOOe+00 



Figure 7: Typical output produced by running afga.m. 

We've written 3 Octave/Matlab m-files called af ga . m, af ga_step . m and af ga_rot . m 
that implement some of the results of this paper. The main file afga.m calls functions 
in afga_step.m and afga_rot.m. 

The first 3 lines of afga.m instantiate the 3 input parameters gO_degs (=7 in 
degrees), del_lam_degs(= AA in degrees), and num_steps (= maximum value of j 
that will be considered, j will range from to num_steps) 

Each time afga.m runs successfully, it outputs a text file called afga.txt. 
FigJT] illustrates a typical afga.txt file. The first 3 lines record the inputs. The next 
line labels the columns of the file. The column labels are 

• j = J 

• gam_j (degs) 

• alp_j (degs) 

• vr_x = [fj]^ 

• vr_y = [fj]y 

• vr_z = [fj]^ 

• VS_X = [Sj]^ 

14 



= 7j in degrees 
= aj in degrees 




Figure 8: Values of 7j and aj obtained with afga.m for 7 = 169.15° and various 
values of A A 



• vs_y = [sj]y 

• VS_Z = [Sj]z 

Following the line of column labels are nuin_steps+l lines of output data. 

In afga.txt, data in each line is separated by a tab. Thus, the full afga.txt 
file can be cut-and-pasted into an Excel spreadsheet or other plotting software in 
order to plot it. 
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■a 



gamiTia(degs)=2 1.15 
del_lam(degs)=180 




DO 
T3 



gammaj 
alpha_J 



gainnia(degs)=2 1.15 
del_lam(degs)=90 



&0 

-a 



gamma_j 
alpha_J 



gamma(degs)=21.15 
del_lam(degs)=45 




gamnia(degs)=2 1.15 
del_lam(degs )=5 




J 



Figure 9: Values of 7j and aj obtained with afga.m for 7 
of AA 



J 

21.15'^ and various values 



FiglH] shows the values of ■jj and aj obtained with afga.m for 7 = 169.15° and 
various values of AA. FigJH shows the same thing but for 7 = 21.15°. We see that 
decreases almost linearly from 7 to near zero. The behavior of near zero depends 
on the value of AA. For < AA < vr, 7^- goes to zero without too many oscillations. 
For AA precisely equal to it, 7^ never converges to zero. It gets trapped near zero, 
oscillating about it with a constant amplitude. 

In Appendix \^ we discuss in more detail the behavior of our AFGA when 
AA = TT. This case most closely resembles the original Grover's algorithm. 

In Appendix [Bl we discuss the continuum limit where A'-fj tends to zero for 
all j. This limit is a smoothed out version of the discrete case. It is easily solved, and 
gives a good idea of the rate of convergence of 7^ (and of ERR) towards zero (when 
AA 7^ tt) as j tends to infinity. 
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A Appendix: When AA = tt 



^l<^ \ 

A) 


/ A 

/ X^* 




fo^^" — »J 

\ \ ^^^/^ 
\ \ / * 

,2 \ Si 



A 

-z 




(a) 



(b) 



Figure 10: (a)sj vectors for original Grover's algorithm. (b)sj vectors for our AFGA 
with A A = TT and the same 7 as in (a). 





(a) 



(b) 



Figure 11: (a)r when 7^,^, > ^. (b)r when 7^-^^, < ^ 



A7 



In this section, we discuss the AA = vr case of our AFGA. 

FigJTUlfa) shows the pattern of the Sj vectors for the original Grover's algo- 
rithm, and Fig JTOT b) shows the pattern for our AFGA with AA = vr and the same 
7 as in (a). There is a critical j, call it jsat ("sat" for saturation), (a) and (b) have 
the same Sj vectors for j = 0, 1, . . . , jsat- For j > jsat-, the Sj vectors of (a) continue 
to decrease their angle (with respect to z) at a uniform rate, past the North Pole, 
whereas the Sj vectors of (b) get trapped in the neighborhood of the North Pole, 
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bouncing back and forth, making an angle of ±r with respect to z. In other words, 
the pattern observed is hke this: 



3 


7i for (a) 


7i for (&) 








jsat - 1 


Ijsat + A7 


7i.at + A7 


3 sat 


Ijsat 




jsat + 1 


Ijsat - A7 


-r 


jsat + 2 


Ijsat - 2A7 


r 


jsat + 3 


lisat - 3A7 


-r 


jsat + 3 


Ijsat - 4A7 


r 









jsaf is defined by the constraints that 7j^„t-i > A7 and < 7^-^^^ < A7. 
As shown by FigJTTl 



7i , if < 7,- < ^ 

ijsat — ijsat — 2 

A7-7,.„. if^<7....<A7 
min(7j,,,,A7-7j-^^J . 



For example. 



7(degs) 


A7(degs) 


7j.at(degs) 


r(degs) 


160 


2(20) 








164 


2(16) 


4 


4 


166 


2(14) 


26 


2 



(92) 
(93) 



(94) 



The last column of Eq. flMj) was calculated using Eq.f l93|) and then checked using 
af ga.m. 



B Appendix: Continuum Limit 

In this section, we explore 7^ in the continuum limit. 

Suppose we take the limit where A7J tends to zero for all j > 0. We replace j 
by a real number t > and ■jj by a continuous function g(t) of t. A'jj = ■jj — 7^+1 — 
and Eq.dZl]) tends to 

- ^ = -7 + ^ + atan2(|5^(g)|,C^(g)) , (95) 

where fi{g) satisfies 

C/x{c/) = C^Cg + S^SgC^x ■ (96) 
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Figure 12: A spherical triangle with sides of length /i, 7 and g, and angle AA between 
the g and 7 sides. The spherical triangle has sides which are segments of great circles 
of the unit sphere. 



Henceforth, we will restrict our attention to the case < (7 < 7 < tt. 

In general, note that if cos(a) = cos(/9), then a = ±/3 + 2ttN for some integer 
A^. Hence, Eg. (1961) doesn't specify fi{g) uniquely. However, the "Law of Cosines" of 
spherical trigonometry tells us that one possible value for fi{g) is the length of the 
side of the spherical triangle portrayed in FigJT2l Henceforth, we will identify fi{g) 
with this unique geometrical value. When /i is given this geometrical value, since 
< < 7 < TT, /i G [0, 27r]. Since atan2(|S'^|, C^) G [0, tt], we must have 

atan2(|S„|.C„) = { _ ^ Zl^M 

= min(/i, 27r — /i) . (98) 
Hence, when ^{g) is given its geometrical value, 

dg 

- = -7 + fl- + mm(/i, 27r- fi) . (99) 
Eq.( p9|) can be solved in closed form in some special cases: 
(a) |5f| << 1: In this case, we get from Eq. (!96!) . 



C,^C, + S.gCAx ~ cos(7 - gCAx) ■ (100) 
so ~ 7 — gC^x- Hence, 

- ^ ^ ^(1 - Caa) (const.)e-*(^-^^^) . (101) 

dt 

(b) < ■y — g « 1: In this case. 



^ ~ min[/i(7),27r-^(7)] ^ ^ ^ 7 - t min[/i(7), 27r - /i(7)] (102) 
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(b.l) AA = vr: In this case, we see from FigJT2] that /i(7) = 27. Hence 

\ =,-tl^-< if o^>| ^ 

^ \ 7 - t27 if 7 < f 

(b.2) < AA « 1: In this case, by virtue of Eq.®, 

C,^C^^ + S^^Ax ~ 1 - ^5'(AA)2 ^ cos(5-,AA) , (104) 
so /i(7) ~ S^AX. Hence 

^ ^ 7 - t^^AA . (105) 
And what is the maximum \dg/dt\ for (7^7? When ~ 7, 

max \dg/dt\ = \dg/dt\Ax=TT = min[27,27r — 27] . (106) 

AAe[0,7r] 
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